First-principles study of lithium-doped carbon clathrates under pressure 
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We present a theoretical study on the behavior under pressure of the two hypothetical C46 and 
LigC46 type-I carbon clathrates in order to bring new informations concerning their synthesis. Using 
ab initio calculations, we have explored the energetic and structural properties under pressure of 
these two carbon based cage-like materials. These low-density meta-stable phases show large nega- 
tive pressure transitions compared to diamond which represent a serious obstacle for their synthesis. 
However, we evidence that a minimum energy barrier can be reached close to 40 GPa, suggesting 
that the synthesis of the Li-clathrate under extreme conditions of pressure and temperature may be 
possible. Electronic band structure with related density of states behavior under pressure as well as 
the dependence of the active Raman modes with pressure are also examined. 
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I. INTRODUCTION 

Hypothetical carbon analogues to silicon clathrates are 
predicted to exhibit promising properties of fundamen- 
tal and technological interest. These low-density class of 
cage-like materials, related to the zeolite topology, were 
first investigated by Nesper et al 1 . They allow the in- 
tercalation of guest atoms inside the cages and therefore 
display unique properties. Exceptional values of mechan- 
ical properties or superconductivity at high temperature 
are predicted within the framework of ab initio calcula- 
tions 1 ' 2,3, . Also the low work function of the metallic 
LigC46 clathrate would be interesting in the design of 
potential electron-emitting devices 5 . Contrarcly to the 
other group-14 clathrates (Si, Ge and Sn), the synthesis 
of the carbon clathrates remains a challenge up to now 
even if the deposition of C20 or C36 networks, which basic 
units have a comparable size to the ones of the carbon 
clathrates, have been prepared 6,7 . 

Type-I semiconducting carbon clathrates consist of the 
assembly face-sharing C20 and C24 nano-cages result- 
ing in a cubic structure (SG: Pm3n, No 223) (Fig. 1). 
They would represent a low-density fourfold diamond- 
like sp 3 configuration. Intercalation with the Li and 
Be small atoms inside each cage of the type-I carbon 
clathrate has been suggested 8 . Ab initio calculations 
indicate that upon efficient doping, carbon clathrates 
would turn metallic 2,5,9,10,11 . Furthermore, large val- 
ues of the electron-phonon interaction potential let en- 
visage carbon clathrates as possible candidates for high- 
temperature superconductivity 4,9,11 . Recently, Zipoli et 
al have calculated a large value of the electron-phonon 
coupling for the clathrate FC34 leading to a supercon- 
ducting temperature close to 77 K 12 . In addition, X. 
Blase et al have shown that the minimum strength of 
the carbon clathrates should be larger than that of di- 
amond 3 . C-clathrates have been also predicted to have 
extraordinary mechanical properties as a bulk modulus 
close to the one of diamond 2,13 or even high strength 
limits under tension and shear beyond diamond ones 3 




FIG. 1: Symbolic representation of the C20 and C24 face- 
sharing carbon nano-cages. 



The silicon typc-I clathrate Ba 8 Si46 can be synthesiscd 
from a precursor consisting of a stoichiometric mixture 
of BaSi 2 and Si at 3 GPa and 800°C 14 . An attempt was 
made by Yamanaka et al to prepare a carbon clathrate 
using high pressure (HP) and high temperature (HT) 
conditions from the polymerization of Cgo but no evi- 
dence of a clathrate-like phase analogous to Si-clathratcs 
were found 15 . Another possible HPHT route has been 
envisaged using a Li-graphite intercalation compound, 
where the Li atoms could frustrate the graphite-diamond 
conversion and promote the formation of the endohc- 
drally doped nano-cages 13,16 . This idea is supported by 
the observation of the half conversion of the sp 2 bonds 
into sp 3 bonds at 17 GPa when graphite is pressurized 17 . 

In this work, we investigate the high-pressure proper- 



ties of the C46 and LisC46 carbon clathrates within the 
framework of the Density Functional Theory as a way to 
obtain new informations concerning their synthesis. The 
diamond-clathrates transition is discussed as well as the 
prospect of a reaction path leading to the Li-intercalated 
carbon clathrates. The electronic band structure proper- 
ties of C46 and Li§C46 under pressure are also examined. 
We have calculated the active Raman modes of LigC46 
and studied their pressure dependence which could con- 
tribute to the identification of carbon clathrate crystals 
if synthesized for instance in a diamond anvil cell. 

II. COMPUTATIONAL DETAILS 

Our ab initio calculations were performed with the 
VASP code 18 , using the Projector Augmented Waves 
(PAW) 19 pseudopotentials supplied with the package 20 , 
within the Generalized Gradient Approximation (GGA) 
for the exchange-correlation energy 21 . This choice is 
justified as in the case of sp 3 carbon based materials, 
GGA gives a better evaluation of the C-C distance 22,23 . 
A plane-wave energy cutoff of 650 eV together with a 
(4x4x4) Monkhorst-Pack grid corresponding to 10 k- 
points in the irreducible part of the Brillouin zone were 
used for the C46 and LigC46 compounds. Our k-mesh 
was enough to reach the convergence due to the large 
clathrate cells. Thus, convergence of 2 meV/atoms in 
the difference of total energies and 0.2 GPa in pressure 
was reached. The cell and internal parameters were fully 
relaxed in order to obtain forces over the atoms lower 
than 0.005 eV/A. The resulting energies as a function 
of the volume were fitted using a Murnaghan equation of 
state (EOS) 24 to obtain the equilibrium volume, the bulk 
modulus and its pressure derivative at zero pressure. The 
pressure calculated values are in excellent agreement with 
the ones provided from the fit. The study of the enthalpy 
allowed us to determine the relative stable phase. 

III. RESULTS 
A. Structural properties 

Figure 2 shows the calculated energy-volume (E-V) 
curves of the carbon phases considered, i.e. diamond (C2) 
which is included for comparison and the C46 and LigC46 
clathrates structures. The diamond phase has the lowest 
energy and is the most stable at zero pressure. In Ta- 
ble I, the results for the three studied structures are given 
from our calculations together with existing experimen- 
tal values. The internal parameters of the equilibrium 
volume of C46 and LigC46 are shown in Table II. Our 
calculated values for the lattice parameter and the bulk 
modulus of diamond are in excellent agreement with the 
experimental works 25,26 . The resulting equilibrium lat- 
tice parameter and bulk modulus for C46 (a=6.696 A, 
.Bo =371 GPa) are close to those from previous LDA and 
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FIG. 2: Calculated Energy- volume curves for diamond (C2), 
C46 and LigC46- Vo la is the atomic volume of the diamond 
structure at equilibrium. 



GGA studies 5,12 . The intercalation of Li atoms at the 
center of the C20 and C24 cages leads to an expansion 
of the lattice parameter (a=6.833 A) as well as a higher 
compressibility (Bo=356 GPa). Atomic volume has been 
slightly expanded by 14.5 % and 3.5 % for the C46 and 
the LigC46 phases respectively with respect to diamond. 
We find an energy difference of 92 meV/atom for C46 
relative to diamond in good agreement with Nesper et 
al (90 meV/atom) 1 who used also the PAW method. 
Other ab initio methods give much higher energy dif- 
ferences 144-210 meV/atom 27,28,29,30 . Despite the dis- 
crepancies found from all the mentioned calculations, it 
is remarkable that these binding energies, relative to di- 
amond, are considerably lower than that of the fullcrcne 
Ceo (430 meV/atom) 31 . However, according to our cal- 
culations, filling the voids inside the C46 clathrate dra- 
matically increases the enthalpy at zero pressure to 1.48 
eV/atom. 

The compressibility of the C46 and LigC46 clathrates 
is close to the diamond one until ~8 GPa (Fig. 3). Then, 
C46 and Li 8 C46 become slightly more compressible than 
diamond, the pure clathrate being less compressible than 
the filled clathrate. This contrasts with the observations 
in Si-clathrates 32 in which intercalation trends to shift 
the bulk modulus towards the Si-diamond phase. The 
pressure-evolution of the C-C and Li-C distances in the 
cages of the two clathrates are displayed in Fig. 4. The 
insertion of lithium inside each cages of C46 leads to a 
slight extension of the cages and produces a general stiff- 
ening in the pressure-evolution of the C-C distances. The 
Li-C and C-C distances decrease smoothly with pressure. 
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TABLE I: Structural and cohesive parameters of diamond 
(C2), C46 and LigC46- The equilibrium volume Vo and cell 
parameter ao are obtained after relaxation of the forces. Bulk 
modulus Bq and its pressure derivative B' Q are obtained with 
the Murnaghan EOS fit. The relative atomic volume V/Vo 
and enthalpy AHo with respect to the diamond phase are 
also given. 





ao 


Vo 


V/Vo 


Bo 


B' 






(A) 


(A 3 ) 




(GPa) 




(eV/atom) 


C-2 


3.576 


11.4 


1.000 


437 


3.22 







3.567° 






442 s 






C46 


6.696 


300.2 


1.145 


371 


3.44 


0.09 


LigC46 


6.833 


319.0 


1.035 


356 


3.13 


1.48 



"Experimental value from reference 26 . 
'Experimental value from reference 25 . 



TABLE II: Internal parameters of the C46 and LigC46 
clathrates at their equilibrium volume. 



C 46 , SG: Pm3n, No 223, Z=l 


Atoms 
C(l) 
C(2) 
C(3) 


Wyckoff positions 
6c 

16i (x, x, x) 
24fc (0, y, z) 


Internal parameters (x, y, z) 
0.25 0.5 
0.1845 0.1845 0.1845 
0.3056 0.1187 


Li 8 C 4 6, SG: Pm3n, No 223, Z=l 


C(l) 
C(2) 
C(3) 
Li (1) 
Li (2) 


6c 

16i (x, x, x) 
24fc (0, y, z) 

2a 

6d 


0.25 0.5 
0.1854 0.1854 0.1854 
0.0345 0.1213 

0.25 0.5 
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FIG. 3: Volumic variation as a function of pressure for (C2), 
C46 and LigC46- 



FIG. 4: Top: Interatomic distances evolution with pressure 
for LigC46 (a): Li-C (top panel) and C-C distances (bottom 
panel), and for C46 (b): C-C distances. 



There is no sudden drop corresponding to a cell volume 
reduction like it has been experimentally observed in the 
Ba-Si and Si-Si distances of Ba 8 Si4 6 33 ' 34 . 



B. Stability under pressure 

The diamond structure is the most stable phase stud- 
ied in this work for positive pressures (Fig. 5). Enthalpy 
curves of C46 and LisC^ are found to be lower than the 
one of diamond at -19 GPa and -187 GPa, respectively. 
Our C2— >CiQ pressure transition agrees reasonably with 
the one calculated by Perottoni and Da Jornada within 
an Hartree-Fock approximation 30 . Negative pressure val- 
ues are expected in those kind of transitions implying ex- 
panded phases compared to diamond. Although, in abso- 
lute values, they are much larger than those calculated in 
the case of the Si 2 ^Si 46 (-6 GPa) 30 and Gc 2 ^Ge 46 (-2.4 
GPa) 35 transitions. In the case of the intercalated type-I 
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FIG. 5: Enthalpies of the diamond (C2) and the C46 and 
LigC46 clathrates phases as a function of pressure. The insets 
provide a zoom of the regions of interest. 
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FIG. 6: Relative enthalpy of LigC.46 from Eq. 2 as a function 
of pressure. AH (LigC^e) >0 for the pressure range studied. 
A minimum corresponding to an energy barrier of 1.9 eV is 
found at 38-39 GPa. 



studied the variation of enthalpy of the previous reaction, 
defined as follows 

AH(Li s C 46 ) = H{Li 8 C 46 ) - 8H(Li) - 23ff (C 2 ). (2) 

Figure 6 shows the variation of AH (LisC 4 e) as a func- 
tion of the pressure. AH(LisC 4 e) is always positive in- 
dicating that the considered reaction 1 is not favourable 
in the pressure range studied. However, we have found 
a minimum corresponding to an energy barrier of 1.9 eV 
around 38-39 GPa. 



C. Electronic properties under pressure 

We have calculated the electronic band structure of 
pure and doped carbon clathrates choosing the same path 
along the same high-symmetry directions of the Brillouin 
zone as in Rcf. 36 were chosen. The band structure of 
the semiconducting C46 at zero pressure and at 51.3 GPa 
is provided in Fig. 7(c) and (d). We find a quasi-direct 
band gap of approximatively 3.88 eV at zero pressure 
between F-X which is in good agreement with previously 
reported results 5,37 . The band gap value decreases with 
pressure, reaching 3.69 eV at 51.3 GPa, and becomes di- 
rect from 2.8 GPa. We can check in Figure 7(a) that 
the intercalation of lithium in each cage of C46 leads to 
a metallic behavior with the filling of the conduction- 
band by Li 2s electrons. A maximum in the density 
of states at the Fermi level Ep is observed suggesting 
interesting superconducting properties as previously re- 
marked 2,4,5,11 . No dramatic changes have been found at 
55 GPa in the band structure of Li§C46, especially in the 
conduction band at the vicinity of Ep (Fig. 7.(d)). In 
this regard, a band crossing over the Fermi level in the 
conduction-band of Ba§Si46 near the R-point, ascribed 
to an electronic topological transition, has been calcu- 
lated recently 34 . In LisC46, the Fermi energy is found to 
increase with pressure, corresponding to a shift of AEp 
~ 1.7 eV between and 55 GPa. This indicates that 
the pressure is a valuable parameter to tune the charge 
transfer in these guest-host systems. We point out that 
the density of states around present a maximum at 
55 GPa. 



D. Dynamical properties of LisC46 under pressure 



carbon clathrates with Na (Pt=-77 GPa) 30 or Li, huge 
negative pressure transitions relative to diamond seem to 
be systematic and render difficult the observation of the 
doped clathrate phases. 

In order to hedge the phase stability of LigC46, we 
considered the following reaction: 

8Li + 23C 2 ^Li 8 C 46 , (1) 

involving stoichiometric ratio of metallic lithium (Li) and 
diamond (C2) to obtain the LisC46 compound. Thus, we 



We used the direct approach based on the small dis- 
placements to calculate the eigenvalues of the dynamical 
matrix at the r-point of the Brillouin zone ' . In order 
to achieve accurate results, we relaxed the cell and the 
internal parameters to obtain forces over the atoms below 
0.0001 cV/A. We performed our study with the PHON 39 
and PHONON 38 programs using a (lxlxl) supercell. In 
order to obtain well converged results, we tested different 
small displacements, and (4x4x4) k-mesh was used to 
reach convergence of the phonon frequencies lower than 
3-4 cm -1 . 
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FIG. 7: Band structure and density of states associated for L18C46 at GPa (a) and at 55 GPa (b) and for C46 at GPa (c) 
and at 51.3 GPa (d).In the case of the LisC46 metallic compound, the Fermi energy is taken as reference. In the semiconducting 
C46, the top of the valence band is chosen as the energy reference. 



TABLE III: Raman active modes of L18C46 (point group: O^) 
for the equilibrium volume (at the T-point). The frequency 
at 870 cm -1 consists in the combination of a Raman (R) and 
an infra-red (IR) mode. These frequencies and their repre- 
sentation were calculated with the phonon program 38 . 

Representation Frequency (cm -1 ) Activity 



T 29 


247 


R 


E 9 


310 


R 


E 9 


387 


R 


T 2 9 


425 


R 


E 9 


550 


R 


T 29 


649 


R 


E 9 


654 


R 


Ai s 


674 


R 


E 9 


690 


R 


T 2 9 


766 


R 


E 9 


804 


R 


T 29 


813 


R 


T 29 


864 


R 


Ai 9 +Ti u 


870 


R + IR 


T 2 9 


921 


R 


E 9 


944 


R 


T 29 


959 


R 


E 9 


976 


R 


T 2 9 


1049 


R 


Alg 


1080 


R 



In order to link our theoretical results with possible fu- 
ture Raman spectroscopy measurements, if synthesised, 
we present the list of the active Raman modes of the 
Li§C46 clathratc corresponding to the equilibrium vol- 
ume in Table III, that could be helpful as a reference for 
the experimental study. 

Figure 8 shows the evolution with pressure of the cal- 
culated Raman frequencies obtained with the phon pro- 
gram 39 . We point out that a typical discrepancy of 
15 cm -1 between the eigenvalues calculated from the two 
programs phon and phonon was observed. This value 
can be considered as an average incertitude of our calcu- 
lated frequencies. The pressure dependence of some fre- 
quencies can be reproduced with a second-order polyno- 
mial fit as shown with the continuous lines in the figure. 
For the other frequencies, a more complicated behavior is 
observed. In that case, we can distinguish three regimes 
with different effective slopes. 



IV. SUMMARY AND CONCLUSIONS 

In this work, we have presented an ab initio theoreti- 
cal study of the behavior under pressure of the C46 and 
Li 8 C46 carbon clathrates. We have analysed the phase 
stability of these two compounds with respect to the dia- 
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TABLE IV: Coefficients of LisC46 Raman frequencies pressure dependence shown in Fig. 8. A second-order polynomial function 
was used to fit the calculated frequencies. 

v(P) = a + bP + cP 2 

Eg Eg T2g 

a (cm -1 ) 263.28 322.83 686.84 1067.68 

b (cm _1 /GPa) 1.69 3.18 3.68 3.87 

c (cm^/CPa 2 ) 1.39 x 10" 2 1.95 x 10 -2 2.26 x 10~ 2 1.45 x 10 -2 




10 20 30 40 50 
Pressure (GPa) 



60 



milestone for future experiments. First, our calculations 
confirmed the interesting structural properties of these 
clathrates. In particular, these structures are of low com- 
pressibility but still more compressible than diamond. As 
a result, no cell volume collapse with increasing pressure 
was evidenced like it occurs in silicon clathrate in the 
same studied pressure range. Our calculations on the 
C46 clathrate are in a good agreement with the previous 
studies. In the case of LisC^, we find a large enthalpy 
at zero pressure related to diamond. Furthermore, the 
Li-intercalated clathrate becomes the stable phase rel- 
ative diamond near -187 GPa, pointing to a difficulty 
to synthesize this compound. We have also proposed of 
reaction path involving stoichiometric ratio of diamond 
and lithium to obtain the LisC46. Although this reaction 
seems not be favourable, we have calculated a minimum 
corresponding to an energy barrier of approximatively 
1.9 eV around 38-39 GPa, an easily accessible pressure- 
window with the diamond anvil cell. Hence, even if the 
phase stability of LisC46 appears to be a serious obsta- 
cle, our result lets suppose that HPHT conditions are well 
suited to obtain this compound. Ab initio cannot pro- 
vide a synthesis route, only future attempts will allow 
to find it. Further ab initio calculations on intercalated 
clathrates considering others chemical species or with dif- 
ferent structures will be valuable in order to reduce the 
cohesive energy. 



FIG. 8: Pressure-evolution of the calculated LisC46 Raman 
frequencies with the phon program 39 . The continuous lines 
represent a polynomial fit (see Table IV). 

mond phase and their structural and electronic properties 
under pressure. We have also calculated the active Ra- 
man frequencies as well as their pressure dependence as a 
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